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We investigate the morphology of a toroidal fluid membrane vesicle confined inside a spherical 
container. The equilibrium shapes are assembled in a geometrical phase diagram as a function of 
scaled area and reduced volume of the membrane. For small area the vesicle can adopt its free form. 

When increasing the area, the membrane cannot avoid contact and touches the confining sphere 
along a circular contact line, which extends to a zone of contact for higher area. The elastic energies 
of the equilibrium shapes are compared to those of their confined counterparts of spherical topology 
to predict under which conditions a topology change is favored energetically. 

PACS numbers: 87.16.D-, 87.10.Pq 


I. INTRODUCTION 


Topological transitions abound in physical and biological systems. In fluid dynamics they can be induced by 
capillary instabilities like the Plateau—Rayleigh instability which causes a fluid thread to break up into smaller 
droplets [T]. Another classic example concerns area-minimizing surfaces: a catenoidal soap film spanning two circular 
wires undergoes a transition to two soap films covering both wires separately when the distance between the wires 
is increased [2]. A topological transition can also be induced by deforming the wire such as in the case of a Mobius 
strip which can change its topology with a twist singularity [3] . 

In the biological world, cells are controlling the topology of their ingredients at all length scales. At the molecular 
level the function of DNA [4] and proteins [5] depends crucially on their topology. On a larger scale membrane-bound 
organelles such as the endoplasmic reticulum or mitochondria display intricate geometries whose topological features 
are carefully tuned for proper function ISH- A healthy cell can change the shape and topology of its membranes 
with the help of other constituents such as the cytoskeleton unmi]. However, topology changes can also occur 
due to pathological conditions: one example is X-linked centronuclear myopathy—a congenital muscle weakness—in 
which the membrane stacks of the sarcoplasmic reticulum of skeletal muscle cells are remodeled into cubic membrane 
structures of higher genus [12] . Similar structures have been observed in the inner membrane of mitochondria upon 
starvation m- 

A theoretical ansatz to model the morphology of the mitochondrium consists of studying a fluid membrane vesicle 
inside a confining spherical container [MlH^ . Although this model is too simple to capture all the details of the 
complex folding patterns, one obtains the basic morphological building blocks. Experimental observations in vivo 
m and in artificial systems m confirm the theoretical predictions. However, the spherical topology of the confined 
membrane vesicle has been fixed in all existing theoretical studies. 

To handle topology changes as well, we will consider confined fluid membrane vesicles of toroidal topology. The 
equilibrium shapes of free toroidal vesicles without confinement have been studied since the 1990s [2TH28] . Only 
recently new stable non-axisymmetric shapes have been found [29j . As in these articles we use the classical curvature 
model to study our system (see Sec. 0 - To obtain a geometrical phase diagram with all equilibrium shapes, we start 
by analyzing the free solutions and determine the parameter values at which the vesicle starts to touch the confining 
sphere (Sec. HI). In the next step we determine the axisymmetric (Sec. IV) and non-axisymmetric (Sec.|V|) equilibrium 
shapes in contact with the confinement. In Sec. |V^we finally compare the vesicles of spherical and toroidal topology 
to predict under which conditions a topological transition can take place. 


II. MODEL 

We search for the equilibrium shapes of a toroidal membrane vesicle of fixed area A and volume V in spherical 
confinement. In the following area and volume are scaled by the respective quantities Aq and Vq of the confining 
sphere: a = A/Aq and v = V/Vb- The elasticity of the membrane is modeled within the scope of the classical 
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FIG. 1: (Color online) Transition from a spherical to a toroidal topology (see text): (bottom left) axisymmetric equilibrium 
solution with spherical topology (a,v) = (1.2, 0.8), (top left) ellipsoid-like equilibrium solution with spherical topology (a,v) = 
(1.5, 0.8), (top right) confined toroidal membrane (a,v) = (1.2, 0.8) (simulation snapshot), and (bottom right) axisymmetric 
equilibrium solution with toroidal topology (a,v) = (1.2, 0.8). The maximum values of the scale bar are Cmax = 10, 23, 11, and 
8 (clockwise starting from bottom left). 


curvature model [3QH32] with the bending energy 

Eh := E, + Eji = '^ j da {2H - Cof + R j dA K , (1) 

where H is the mean curvature and K the Gaussian curvature. The stiffness of the membrane is reflected in the 
material parameters k and R which are the bending and the Gaussian rigidity, respectively. We assume in particular 
that the membrane prefers to be intrinsically flat which implies that the spontaneous curvature Cq is zero in our 
study. The second term of the energy involving the Gaussian curvature is a topological invariant for a closed surface. 
According to the Gauss-Bonnet theorem it is given by [33] 

Ek = R J dA K = 47r^(l — g) (2) 

for a surface of genus g, where g corresponds to the number of “handles” of the vesicle. For a spherical surface such as 
those studied in [ans] g is zero; for the toroidal surfaces considered here g is equal to one and Ej^ vanishes. Depending 
on the sign of R an increase of g will thus correspond to an abrupt increase or decrease in energy, respectively. For 
fluid lipid mono- and bilayers one typically finds R ^ —k < 0 [34], which implies that a change from a spherical to a 
toroidal topology will increase the term of the energy involving K. However, the first term of the energy will change 
as well and it is not clear a priori whether the total energy will increase or decrease. 

Why should the topology of the membrane change at all? The equilibrium shapes of a spherical membrane vesicle 
in confinement have been determined recently by some of the authors nuns]. For moderate values of area (a > 1) 
we have found an axisymmetric solution with a “light bulb”-like invagination which starts to contact itself when the 
area of the vesicle is increased (see left part of Fig. [^. If the membrane was able to fuse, for instance, with the 
help of embedded specialized molecules such as SNARE proteins |35|, the two parts of the surface in contact could 
merge. The topology of the resulting surface is toroidal (see right part of Fig. [^: two necks instead of one connect 
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FIG. 2: (Color online) Partial phase diagram {a^Vr) showing the free axisymmetric equilibrium shapes together with the 
conformal transformations of the Clifford torus. The colored regions indicate the global minima of this set of shapes, which are 
either sickle-shaped tori (lime green region), circular tori (orange region), or conformal transformations of the Clifford torus 
(hatched orange region). The discoid tori (purple) are local minima of higher energy. The curves ^ delimit the regions 
in which the free shapes are equilibrium solutions of the confinement problem. To their right no such solution exists, which 
explains why parts of the phase diagram are white. Inset: Scaled bending energy := Eh/(STiHi) as a function of Vr for the 
different solution types (see Ref. [23l[27]): sickle-shaped (lime green curve ^), circular (red curve [2]), discoid (purple curve 
[4]), and stomatocyte (dark green curve [^) tori, and the conformal transformations of the Clifford torus (blue curve ^). The 
long-dashed parts of these curves correspond to shapes of maximum bending energy. The short-dashed line ^ corresponds to 
the limit torus with = 2, which consists of two identical spheres connected via two infinitesimal holes. 


the invagation to the outer part of the membrane in contact with the spherical container. The toroidal ground state 
for the parameters of Fig. is axisymmetric. We will see in the following that the symmetry of the solution depends 
on the respective values of area and volume of the vesicle. 


III. FREE SOLUTIONS 

Because of the scale invariance of the bending energy the only relevant parameters of the free vesicle are the 
reduced volume Vr = and the spontaneous curvature Cq. For vanishing spontaneous curvature the following 
types of axisymmetric local extrema have been found (see Refs. [umiiiiT] and Fig. [^: (i) sickle-shaped tori, {ii) 
(almost) circular tori, {iii) discoid tori, and {iv) toroidal stomatocytes. All but the stomatocytes have an additional 
reflection symmetry with respect to the plane perpendicular to the axis of symmetry. 

The lower limits of discoid tori and toroidal stomatocytes are shapes with a vanishing hole diameter at the vertical 
axis (purple and dark green squares in Fig.[^. The lower limit of the stable sickle-shaped tori, 0, corresponds to 
two connected identical spheres with an infinitesimal hole in each of the poles. In the same limit the (almost) circular 
tori approach a shape of exact circular cross section while both the diameter of the hole and the cross section vanish 
[23]. At Vr = = 25 / 4 ^ 1/2 ~ 0.71 the cross section of the corresponding profile of the same branch is exactly circular 

as well (blue circles in Fig.|^. This is the Clifford torus whose bending energy is a lower bound to that of any toroidal 
shape [36] . Its two different ellipsoidal deformations have been identified as discoid Vr < vf and sickie-shaped Vr > vf 
tori m- In this articie we wiii use a siightiy different definition by caiiing only those tori sickle-shaped, which have a 
sickle-shaped cross section, i.e., a cross-sectional meridional curvature which changes sign twice. With this definition 
the branch of circuiar tori becomes sickie-shaped for Vr > 0.79 (red circies in Fig.|^. Ali other soiutions of this branch 
(and the corresponding soiutions in partiai or fuii contact with the container, see beiow) wiii be caiied circuiar tori in 
the foiiowing. 

In the inset of Fig. [^the scaied bending energy := Eiy/{^Tin) is shown as a function of Vr [22]. The branch 
of the toroidai stomatocytes contains no minima and wiii thus be omitted in the foiiowing. Without the confinement 
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constraint the global minima can be directly read off: for Vr < 0.57 the tori are sickle-shaped. Above this value the 
circular tori have a lower bending energy. Increasing Vr leads to a discontinuous shape transformation (red squares 
in Fig. since the branch of locally stable sickle-shaped tori runs up to = 0.65 (lime green squares). For Vr > vf 
the ground states are non-axisymmetric shapes. They can be constructed with the help of conformal transformations 


Special conformal transformations of the form 


X' = 


(^+A)2 ■ 


( 3 ) 


conserve the toroidal topology of the surface G and leave the bending energy invariant. They consist 

of an inversion about the origin, followed by a translation along a constant vector A, and a second inversion about 
the origin. The resulting surface has a reduced volume which is larger or equal to the reduced volume 

of the original surface. X' can always be rescaled to stay inside the container. Above the global minima are 
thus non-axisymmetric conformal transformations of the Clifford torus since their bending energy = 7r/2 (blue 
line ^ in the inset of Fig. ^ is always lower than that of the axisymmetric solution of same reduced volume (red 
curve in the inset of Fig.^ [HlllIlT]. The most general special conformal transformation ^ involves the three 
parameters A = (A^,, Ay^Az)7 For an axisymmetric surface X only two parameters remain: without loss of generality 
we can orient A to he in the plane defined by the axis of symmetry 2 : and the perpendicular axis x. Since conformal 
transformations of the Clifford torus along the axis of symmetry rescale the shape, there is only a one-parameter 
family of special conformal transformations leading to the global non-axisymmetric minima discussed above m- 

As long as the area of the vesicle is small enough, the confined vesicle does not touch the container at all and 
adopts the shape of the free system. To predict the behavior for larger area, one can rescale the free axisymmetric 
solutions up to the point at which they come into contact with the spherical container. These solutions have been 
obtained numerically with a shooting method (see appendix for the details of the method). They correspond to 
the curves ^ ^ delimiting the parts in the phase diagram^, 1 ;^) in which the free vesicles are the local minima 

(see Fig. [^. To understand how the global energy minima (colored regions of Fig. are found, consider, for example 
Vr = 0.6. Without the confinement the free circular torus is the global minimum whereas the free sickle-shaped 
torus is a local minimum. At a = 0.6 both solutions can be put inside the container without any contact between 
membrane and confinement. However, when we increase the parameter a keeping Vr constant such that we cross the 
red line in Fig. we do not find any circular tori which are not in contact with the confinement. The only free 
solution is the corresponding sickle-shaped torus which is thus the global and only minimum of all free shapes. In 
general, one observes that the sickle-shaped tori start contacting the confinement for an area which is much larger 
than the corresponding value for the circular tori. Consequently, the sickle-shaped tori are the only free minima for 
0 < 'T’r < 0-65 provided the area is large enough. In contrast to that we can neglect the locally stable discoid solutions 
for the discussion in the following since they start contacting the container for even smaller area than the circular 
tori. 

To find the delimiting curve of the region of the free non-axisymmetric solutions, i.e., the blue curve ^ of the phase 
diagram in Fig. we consider the Clifford torus in contact in the horizontal (xy) plane. Its defining radii are given 
by r = - 7 =—r and R = \/2r. In polar coordinates {0, (j)) its shape can be written as: 


x{eA) 


\/2 -|- 1 


((^/2 + sin 0) cos 0, (a/ 2 + sin 0) sin 0, cos 0) 


( 4 ) 


on which we perform the conformal transformation, Eq. With an appropriate rescaling and a final translation 
along A 


X = X'(l - A^) + A , 


( 5 ) 


every point of the inner membrane which lies on the confining sphere will stay on the sphere. Choosing A = AxX one 
finds a surface with the same circular contact line as the original Clifford torus. In the limit A^ 1, {a^Vr) (1,1) 
and the conformally transformed torus becomes a sphere with an infinitesimal handle. The resulting shapes are Dupin 
cyclides whose centers of curvature are located on two confocal conics: an ellipse and a hyperbola ISSIET]. For the 
Clifford torus the ellipse corresponds to a circle of radius R and the hyperbola becomes a line which coincides with 
the axis of symmetry. Using the results of Refs. [Ml EH one obtains for the area and volume of the transformed torus: 


a = 4(1 — A)A[f (e^)-h 


- 1 


= SA-" 


A2 




2 
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(l + 5d^)K(e“)| 
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FIG. 3: (Color online) Equilibrium profiles for fixed reduced volume Vr — 0.7 and increasing area (from left to right) a = 0.5, 
0.8, 0.9, 1.1, and 1.15. The background colors and labels are in accordance with the final phase diagram (Fig. [^. The color 
code of the profiles is the same as in Fig. with emax = 10, 5.5, 6, 10, and 10.5. 

where /C(e^) and £{e^) are the complete elliptic integrals of the first and second kind, respectively [38]. The parameter 
g _ ^ {0,1} is the eccentricity and A = ^ {0.5, i7} is the semi-major axis of the 

generating ellipse. 

Fig-i does not show the whole story yet. We have not yet discussed the conformal transformations of the sickle¬ 
shaped tori. The resulting surfaces have to be compared energetically to the solutions which are in contact with the 
confinement. We will do this in two steps. We first identify all axisymmetric solutions of the confinement problem. 
In a second step we determine the non-axisymmetric confined shapes and assemble everything in one phase diagram. 


IV. CONFINED AXISYMMETRIC SOLUTIONS 


When you press a cylindrical sheet of paper gently against a wall, the contact between the two will be a line. 
The curvature of the sheet perpendicular to the line of contact will generally not be zero. Increasing the force will 
decrease this curvature until it is equal to zero, i.e., equal to the curvature of the flat wall. Above this threshold a 
contact area will start to form. The same phenomenon is observed in our system (see the first four profiles in Fig. 
for an example): increasing the area for a fixed reduced volume one reaches a point at which the free vesicle starts 
to touch the spherical confinement. The corresponding curves in the phase diagram {a,Vr) have been discussed in 
the previous section. Interestingly, all these solutions are in contact with the container in one circle which lies in 
the horizontal symmetry plane. Their curvature perpendicular to the contact line, Kj_^ is always larger than that of 
the sphere, Kj_ = 1 . Increasing the area further thus leads to solutions which are still confined by the circle (called 
contact-circle solutions in the following). Above a certain threshold area, and the circular contact line 

becomes a spherical segment leading to contact-area solutions. 

The corresponding confined axisymmetric solutions have been determined numerically with the same shooting 
method as in the previous section (see appendix]^. The appropriate boundary conditions depend on the case under 
consideration and are summarized in Tab. [TT| of the appendix. Fig. shows the resulting phase diagram together 
with a comparison of bending energies. The dark green line Ci in the diagram separates the region of contact-circle 
solutions from the region of contact-area solutions 
tori 
dashes 


For Vr < 0.65 these solutions coexist with the free sickle-shaped 
Comparing the respective bending energies one finds a discontinous shape transition (lime green line with long 
^^ in the phase diagram), below which the free sickle-shaped tori are the global minima. A first estimate of 
the upper boundary of the region of confined axisymmetric tori (short-dashed lime green curve ) can be found by 
comparing the bending energies of these solutions with those of the conformal transformations of the free tori. We 
would like to stress that these free shapes are not equilibrium solutions in general but serve as a means to determine the 
regions of the phase diagram in which non-axisymmetric shapes with a lower bending energy than the axisymmetric 
solutions exist. We will see below that the curve including all solutions with and without axisymmetry (black curve 
D 2 in Fig. m see next section) lies below the boundary discussed here. 

Consider nrst the free sickle-shaped torus with a reduced volume of Vr = 0.65 (lime green square). It has the 
smallest bending energy of the whole branch ^ of sickle-shaped tori (see Fig. 4(b)). A conformal transformation of 


this shape along the x axis will not change the energy but increase the reduced volume (black dashed-dotted line [^). 
With a rescaling of such a shape or a subsequent conformal transformation along z (which both do not change the 
energy as well) every point of the phase diagram left of the curve can be reached. The resulting surface has to 
be compared to the other equilibrium shapes of same {a^Vr). Along the black dashed-dotted line one finds that the 
energy of the conformally transformed free sickle-shaped torus is lower than the energy of the axisymmetric contact- 
area solution for Vr > 0.76 (black square). Repeating this comparison for smaller area, one obtains the short-dashed 





















6 



FIG. 4: (Color online) (a) Partial phase diagram (a,Vr) showing the axisymmetric extrema together with the conformal 
transformations of the free vesicles. The colored regions indicate the global minima of this set of shapes. Regions of sickle¬ 
shaped tori are shown in green, regions of circular tori in orange/red. Non-axisymmetric regions are hatched. The colors vary 
from light to dark with increasing conhnement. The curve |Ci | separates the region of contact-circle solutions from the region of 
contact-area solutions, whereas IC 2 1 is the boundary between the confined circular and sickle-shaped solutions (for a description 
of the curves |Di | , D 2 , and see text), (b) Scaled bending energy as a function of Vr for the different axisymmetric 
solutions. In addition to the free solutions (solid lime green red [2], and blue curves ^) lines of constant area a are shown 
as well (dashed yellow J] (cl = 0.7), purple (a = 0.9), and pink |III| (a = 1.1) curves). The short-dashed line ^ corresponds 
to the limit torus with eb = 2, which consists of two identical spheres connected via two infinitesimal holes. 


lime green curve D 2 along which the energy of the conformally transformed free sickle-shaped torus equals that of 


the corresponding axisymmetric contact-area/contact-line solution. The curve ends in the point at which it intersects 
with the branch [2] of the free solutions (lime green circle). The part of the curve to the right of the black square 
can be found by comparing the conformal transformations of all other free tori of the sickle-shaped branch [jj with 
the axisymmetric contact-area solutions. Both curves, D 2 and end in the point (a, Vr) = (2, 0), which corresponds 


to a surface consisting of two connected ident ical spheres with an infinitesimal hole in each of the poles (see Sec. III). 
The short-dashed black curve ^ in Fig. 4(a) [ indicates the conformal transformations of the limit shape along x and 
delimits the regions of shapes with and without self-contact (see next section). The analytical expression for this 
curve is given by: 


Vr = 


l-(a-l)i 


3 

a2 


( 7 ) 


To understand the region of confined tori in the center of the phase diagram, consider curves of constant area a (curves 
^ - I III I in Fig. [^. For a small reduced volume Vr the solutions are free sickle-shaped tori which can be mapped on 
one an other via a simple rescaling. Their bending energy does not depend on a and corresponds to the curve U 
in Fig. 4(b) When we cross the curve |Di in the phase diagram (yellow, purple, and pink squares, respectively), a 


discontmous shape transition to a confined state is energetically favorable. Depending on the value of a the membrane 
can either adopt a contact-circle (curve ^) or a contact-area (curves ^ and |III| ) solution. In the former case (a = 0.7), 
the solution adopts a free shape again above Vr = 0.63 (yellow circle). The corresponding be nding energy is given by 


the curves (between yellow and blue circle) and ^ (to the right of the blue circle) in Fig. 4(b)[ In the latter case 
(a = 0.9 and a = 1.1), we can observe in more detail how the (preliminary) upper limit D 2 of the region of confined 
axisymmetric solutions has been obtained. The respective bending energies of the confined solutions (curves 2ll and 
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(b) (c) 

FIG. 5: (Color online) Comparison between (a) the numerical Runge-Kutta result (black dotted) and the result of a finite 
element simulation (red solid) for (a,i’r) = (1.1, 0.7), and (b) the conformal transformation of minimal bending energy (black 
dotted) and the result of a hnite element simulation (red solid) for {a^Vr) — (1.15,0.7) (see Fig. for a three-dimensional 
representation of the shapes). The red dotted curve outside of the container indicates position and relative value of the 
container force |/^| in the finite element simulations, (c) Scaled bending energy Ch and container force density /o at the 
contact line as a function of a for fixed reduced volume Vr = 0.7. The green dotted line corresponds to non-axisymmetric 
contact-area solutions and was obtained with the help of the finite element simulations. All other curves were determined 
with the Runge-Kutta method. Colors and labels are in accordance with the final phase diagram (Fig.|^. The black triangles 
correspond to the shapes shown in Fig. 


|III| ) decrease first before increasing up to the point where the energy is equal to one of a conformally transformed 
free solution (purple and pink circles on D 2 , respectively). One would thus expect a transition to a non-axisymmetric 


free torus for higher Vr. In reality, however, the membrane will not become a free torus again but prefers to adopt a 
non-axisymmetric confined shape as we will see when including all solutions in our discussion. 


V. CONFINED NON-AXISYMMETRIC SOLUTIONS 


As a starting point we consider conformal transformations given by the combination of the equations <© and 
and apply them on the confined axisymmetric solutions X of the previous section: 


X = 


X 

|X|2 


X 

|X|2 


■AP 


(1 — A^) -f A , with A = cst G 


( 8 ) 


Every point of the membrane in contact with the spherical confinement is mapped on some point of the sphere again 
via Eq. <©• The resulting contact surface/contact line is different from the original spherical segment/circle. However, 
one can show that the local curvature condition at the boundary/boundaries of the free part of the membrane is still 
fulfilled (see appendix [B|) . To be sure that the obtained solutions are in equilibrium, one also has to check whether 
the total force on the membrane exerted by the container is zero. This is the case for the axisymmetric solutions X, 
where the container exerts a constant normal delta force in each point of the contact line(s) to equilibrate the system. 
The corresponding container force densi ty fn follows from the boundary conditions of the Runge -Kutta scheme and is 
given by /o = '^o (see appendix Eq. (Bll), and—for an example—the lower part of Eig. 5(c)). A closer look at the 
conformal transformations X of the axisymmetric shapes reveals, however, that the total force on the transformed 
membrane exerted by the container is not necessarily zero (see appendix [b| . 

In analogy to Ref. [27] one might, nevertheless, hope that the resulting surfaces are close to the true equilibrium 
non-axisymmetric solutions. Assuming this for the moment we can compare bending energies in the same manner as 




















FIG. 6: (Color online) (a) Phase diagram {a,Vr) showing all global minima. Regions of sickle-shaped tori are shown in green, 
regions of circular tori in orange/red. Non-axisymmetric regions are hatched. The black triangles correspond to the shapes 
shown in Fig. 3 The curve |D 2 | separates the region of conhned axisymmetric and conhned non-axisymmetric solutions. For 
a de script ion or the other curves, see Fig. (b) Scaled bending energy et as a function of Vr. This figure is equivalent to 
Fig. |4(b)] except that we have included the confined non-axisymmetric solutions as well. 


in the previous section. For every confined axisymmetric solution one can identify the region in the phase diagram 
in which its conformal transformations can be found. To identify the shape of minimal bending energy for a fixed 
{a^Vr) in this set of solutions, one thus has to identify all axisymmetric shapes (of different which can be 

transformed in such a way that the new surface has the area and reduced volume {a^Vr). Comparing the energies 
of all these shapes with the axisymmetric solution of same (a^Vr) allows to determine an approximate boundary D 2 
which separates the region of axisymmetric and non-axisymmetric global minimum solutions. 

To check how well the obtained solutions approximate the real ones, we have searched for equilibrium shapes 
with the help of numerical simulations. The respective finite element method has already been used for confined 
membranes of spherical topology [niiis]. In the simulations the membrane is discretized and evolved according to 
Newton’s equations of motion until a balance of forces is reached (see appendix [C|) . To find an equilibrium surface 
for a fixed (a, i;^), we have used a discretization of the obtained shape of minimal energy as initial mesh. It turns 
out that the axisymmetric solutions and the corresponding results of the finite ele ment s imulations coincide within 
the numerical error which is about 10“^ in scaled units for all parameters (see Fig. 5(a)). The delta forces exerted 
by the container are smeared out slightly but situated at the contact lines as expected. The necks diffuse away from 
the axisymmetric ground state very easily for non-vanishing noise during the simulations, similar to what has been 
observed for unconfined toroidal vesicles before [29]. The non-axisymmetric solutions, and consequently the phase 
boundary, are approximated remarkably well by the respective conformal transformations (with an error of about 
10“^ in 65 for shapes close to the phase boundary). In additio n to t he smeared-out delta forces at the contact lines 
one now observes additional forces at the contact area (see Fig. |5(b)| ). 

In summary one obtains Fig.j^which shows all global minima assembled in one figure. We can now observe what has 
been anticipated in the previous section: a membrane of large constant area (curves 2E and |III| ) breaks axisymmetry 
but stays in contact with the container when we cross the curve ID 2 I towards higher Vr. Depending on the value of 


a the membrane either adopts the shape of a free conformal transformation of the Clifford torus when crossing the 
curve (purple triangle) or tends towards the limit torus, which consists of two identical spheres connected via two 
infinitesimal holes (pink triangle on curve [l7|). 

It is instructive to consider the case of constant reduced volume as well. The upper part of Fig. 


5(c) show s th e 


scaled bending energy as a function of membrane area for a reduced volume of Vr = 0.7 (compare Fig. |^in Sec. IV). 
For small area the membrane is not in contact with the confining sphere and the bending energy does not depend on 
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{a,v) 



toroidal membrane shape 

(1.1,0.9) 

1.91 

1.88 

non-axi. confined 

(1.2,0.9) 

1.97 

1.97 

non-axi. confined 

(1.1,0.8) 

1.95 

1.87 

axi. confined 

(1.2,0.8) 

1.98 

1.94 

axi. confined 

(1.3,0.8) 

1.99 

1.97 

non-axi. confined 

(1.1,0.7) 

1.99 

1.94 

not confined 

(1.2,0.7) 

1.99 

1.97 

not confined 

(1.3,0.7) 

2.00 

1.98 

not confined 

(1.4,0.7) 

2.00 

1.98 

axi. confined 


TABLE I: Comparison of energies of a vesicle of spherical (from Ref. ^4]) and toroidal topology as a function of scaled area 
a and volume v. 


a. When increasing a contact emerges and the energy starts increasing monotonically. We first find contact-circle, 
then c ontact -area solutions. At a = 1.12 the membrane can break axisymmetry to lower its energy (dotted green line 
in Fig. 5(c)) and adopts a confined non-axisymmetric shape. 

For larger values of fixed reduced volume Vr a similar evolution is observed: for small area a the solution corresponds 
to a conformal transformation of the Clifford torus (hatched orange region). By increasing a we cross the blue curve 
defined by Eqs. (H- The non-axisymmetric solutions in this region are in contact with the container in one horizontal 
circle (hatched light red region). When crossing the left black solid line the contact line extends to a zone of contact. 
The corresponding non-axisymmetric tori are first circular (hatched dark red region) and, upon increasing a even 
further, sickle-shaped (hatched dark green region). All the shapes lying to the right of the curve 0 exhibit self¬ 
contacts and are not studied in the present paper. 


VI. COMPARISON TO SPHERICAL TOPOLOGY 

Knowing the global minima of the confined toroidal membrane vesicles we are now in a position to compare its 
bending energy 0 to that of a confined spherical membrane of same area and volume. For a better comparison we 
again scale the energy with the constant Stt^c, Le., the energy of a perfect spherical vesicle: = E}j/{S7rK) = 

For the free membrane vesicles one finds that the toroidal ground state has a lower energy e,^ than the corresponding 
spherical one as long as 0 < < 0.72 (see Ref. [23]). For larger values of Vr the opposite is true. The second term 

in the energy, e/^, equals for spherical vesicles, whereas it is zero for vesicles of toroidal topology (see Sec. ffl. 
Absolute value and sign of ^ will decide which topology is preferred energetically. Taking a typical value for fund 
lipid mono- and bilayers ^ — 1 [34| the systems increases its energy by Aej^ ^ | when passing from a spherical to 

a toroidal topology. Since this energy is not compensated by the change Ae^ (compare inset of Fig. with Fig. 8 of 
Ref. [39|) we find that the spherical topology is always preferred for the free system. 

A similar comparison can be drawn for the confined vesicles. We focus on the case a > 1 for which we have found 
axisymmetric solutions of spherical topology [HITS] such as the one shown in Fig. [T] {b ottom left). Table depicts 
the values of for both topologies. For a better comparison with the results of Refs. [TH [T5| we have fixed v instead 
of Vr. One observes that the system decreases its energy when changing its topology from spherical to toroidal. 
However, as is the case of the free system, this change does not compensate the increase of ~ ^ in the Gaussian 
term of the energy. Note that the situation changes dramatically if we allow for other values of R. For positive R 
the toroidal topology is favored energetically implying the possibility of a topological transition without any loss of 
energy. 


VII. CONCLUSION 

In this article we have studied the morphology of toroidal fluid membrane vesicles in spherical confinement. A 
combination of analytical theory, numerical calculations and finite element simulations allowed us to determine the 
ground states as a function of area and volume of the vesicle. In contrast to confined membrane vesicles of spherical 
topology we have found non-axisymmetric shapes which do not exhibit any self-contacts. A comparison of energies 
revealed, that the spherical topology is preferred for typical values of the material parameters. 
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FIG. 7: Parametrization of an axisymmetric toroidal membrane vesicle in spherical confinement. The membrane consists of a 
spherical segment in contact with the container and a free part which is given by the solution of the Hamilton equations (A3) 
together with the appropriate boundary conditions (Eqs. d^/Tab.im. 


However, this study was performed assuming zero temperature. To be closer to the real biological world, one would 
have to include the effects of thermal fluctuations. In extensive bilayer fluid membrane phases they can significantly 
affect the competition between spherical morphologies and surfaces of higher genus j4Qf[42] . In these cases an attractive 
interaction appears in addition to Helfrich’s steric repulsion [43]. Guided by experiments on real mitochondria m, 
we expect a similar competition in a biologically significant setting, whenever the membrane is in self-contact. This 
gives a hint to what can happen in the region of self-contacts in the final phase diagram (Fig. 6): the toroidal vesicle 
touches itself, inducing eventually a topology change towards a higher genus (“a sphere with two handles”). 

A study of these shapes and the effects of the temperature goes beyond the scope of this paper but remains an 
interesting open problem for the future. The present model can as well be extended to incorporate a constant adhesion 
energy density due to the contact between membrane and confinement. 
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Appendix A: Numerical method for the axisymmetric case 


The details of the numerical method we have used to find axisymmetric equilibrium solutions can be found in 
appendix A of Ref. m and the references cited therein. The confined membrane consists of two parts, one which is 
in contact with the spherical container of unit size and one which bulges inside to minimize its elastic energy. We 
parametrize the shape of the membrane in terms of the angle-arc length parametrization '^(s) (see Fig. [^, where 
is the angle between the horizontal axis and the tangent of the membrane. The equation describing the evolution of 
^p{s) along the cross section of the free part is an ordinary differential equation of fourth order which can be rewritten 
with the help of Hamilton’s formalism. 

The scaled energy of the free membrane is given by: 


E = E/iiTK) = f dsL 



2 

+ 2dp + Pp^ sin7/^ + Ap(p — cost/^) + \z{z 


sin?/;) , 


(Al) 


where d = ajn and P = P/n are the scaled surface tension and pressure difference, respectively. For the Clifford 
torus one finds d = P = 0. The Lagrange multiplier functions \p and A;^ couple the angle-arc length parametrization 
to the Euclidean coordinate system {p^z). The conjugate momenta of the system are: pp = 2p{ip + 








11 



Section 
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a 

'ipo 
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Pz 

Pp 

P^ 
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IV 


1 

0 

TT 

2 
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55 

55 


55 

Transition contact line ^ area 


IV 


1 

0 

TT 

2 

1 

bisection 

V 

55 



55 

Contact area 


IV 


sin a 

— COSO 

bisection 

1 

see text 

scan 

via Eq. ( 

A4a 

) 

4sino 


TABLE II: Boundary conditions for the different cases discussed in this article. 


Pp = and Pz = Xz- la contrast to the case considered in Ref. [15] the momentum Pz does not vanish 
Hamiltonian is thus: 


i ^ 

yielding the following shape equations 




+ Pp cos 7 /; + sin 7 /; — 2dp — Pp^ sin ip , 


• p^ sm ip 

^ ^ - . 

2p p 

p = cos 7 /^ , 

i = sin ip , 

Pp = ppp‘^ -Pz] cos^p+PpSimp , 


Pp = 




p \Ap 


+ 2d + 2Ppsin'0 , 


= 0 , 

which are solved using a standard shooting method [45] subject to the boundary conditions: 


H 

pU) 

^{s) 


0 , 
Po : 

a . 


^( 5 ) = zo 


]. The 


(A2) 


(A3a) 

(A3b) 

(A3c) 

(A3d) 


(A3e) 

(A3f) 


(A4a) 

(A4b) 

(A4c) 


The values po, zq, and a depend on the case under consideration (see Tab. |^. It is sufficient to scan a and P to 
obtain the free solutions with xy symmetry as well as the transition line of shapes where the contact line extends to 
a zone of contact. For the other cases a third parameter has to be scanned as well. In the case of the free toroidal 
stomatocytes and the contact-circle solutions we have fixed the initial value of ^p. For the c ontac t-area solutions it 
was more practical to scan pz instead. The value of 'ipo can then be obtained by deriving Eq. (A3a), equating it with 
Eq. (A3d) and inserting the values of the other parameters in the resulting equation. 

In all cases the values of 7 /^, p, z, and the corresponding conj ugat e momenta at 5 = 5 are determined by the boundary 
conditions for given a, 'ipQ, and ipQ. The Hamilton equations (A3) can then be integrated with a fourth-order Runge- 


Kutta method. The integration is stopped when the curve returns to the spherical confinement. A bisection on either 
a, ^0 or ^0 (depending on the case, see Tab. O again) determines the correct profiles for which ^{s) = ^ 


— Of. 


Note that we have to fix pressure and surface tension for this method. Area and volume can only be calculated 
for a known profile a posteriori. The equations A = 2iTp and V = 7rp^ simp have to be integrated from s = s to 5 
to find the contribution to area and volume of the free part of the membrane. Adding the result to the contribution 
of the part in contact with the container, Ac = —47rcosa and Vc = |7r(—3cos(a + cos^a), finally yields the values 
for the whole membrane. These allow to assemble the obtained shapes in the geometrical phase diagrams (Figs.[^[^ 
and[^. The Runge-Kutta integration is performed with an integration step-size h = 10“^ in reduced units implying 
a numerical error of 10“^ of the values of a and Vr. 
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Appendix B: Conformal transformations of confined axisymmetric solutions 

We consider the conformal transformation of an axisymmetric surface X given by Eq. (§: 


X 

|X|2 


l|X|2 


•AP 


(1 — A^) + A , with A = cst G 


yielding a deformed surface X which is only axisymmetric when A is parallel to the symmetry axis. The tangent 
vectors of the transformed surface are given by (see chapter 20 of Ref. [46] ) 


Ca = daX = 


|Xo + A|2 


0 2[(Xo + A).eO](Xo+A) 


1 X 0 +A|2 


with Xo = ^ and = daXo = (e^ — 2(X • ea)X)/|Xp. The metric of the transformed surface follows 


Qab — 


(1-A2)2 Q _ (1-A2)2 

|a:o + a| 4 “ |^ + A|4|x 


Qab • 


We orient the normal towards the outside of the vesicle. The transformed normal is then given by: 

2[no-(Xo + A)](Xo + A) 


n = —no 


|Xo+A|2 


(Bl) 


(B2) 


(B3) 


where no = —n + The extrinsic curvature tensor Kab = —n ■ daBb can be obtained by calculating daBb and 

projecting it onto the normal 


Kab = 


|Xo + A|2 


-K 


0 , 2gg,[no-(Xo + A)] ) 

|W + A|2 / 


(B4) 


When the coordinate lines on the original surface before the transformation coincide with the lines of curvature, first 
and second fundamental form are diagonal. These lines will be lines of curvature on the transformed surface as well 
and the principal curvatures follow as {i G {1,2}): 


^ |Xo+A|2 


gu (1 - A2) 
where if? = —\X\'^Ki + 2n-X. The derivative is given by 

1 


0 _L 2[no ■ (Xq + A)] } 
1 X 0 +A|2 / • 


K» + 


dn.Ki 


(1-A2) 


{-|Xo + A|2a„if{> + 2(Xo + A) • (if^o bI - Kf 6°)} , 


(B5) 


(B6) 


where we have used the Weingarten equations in the second step. 

To describe the original axisymmetric surface, we choose the coordinates . The coordinate lines are lines of 

curvature with principal curvatures Ki = K± = ip and K 2 = K\\ = sinT/^/p. For each portion of the surface whi ch is 
in contact with the confining sphere (including the circular contact line(s)) one has = 1. Following Eq. (B5) 

the transformed principal curvatures at the contact lines can be obtained very easily. On the sphere one has |Xp = 1, 
n = X and thus = 1. Moreover, Xq = X, no = n, and one finds that the local contact curvature condition at 
the contact line is also fulfilled for the transformed shape: 


Ki = 


[-|X + A|2 + 2n- (X + A)] = 


1 


-1 - 2X • A - A^ + 2(1 + n • A)] =1 


(1-A2) ... ^ (1-A2) 

The tangent and normal vectors on the contact line are given by b^ = Ba, n^ = n, and 


Ba = 


(1-A^) 
n + Ap 


ea - 


2(A ■ ea){n + A) 
In + Ap 


n = —n + 


2(1 + n • A)(n + A) 

In + Ap 


(B7) 


(B8) 
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The metric is diagonal with gab = |^^|4 9ab- Taking into account that as well as are diagonal, — 

= {K^ — K^)ea = 0. Therefore, the last term in Eq. (B6) vanishes, and one obtains for the derivative of the 
principal curvatures 


daKi = 


- n 


-^Al^daKf |n + A|^ 


(1-A2) 




daKi 


(B9) 


To determine the external forces due to the confinement, consider the stress tensor of the fluid membrane (normalized 
by k) [47] 


r = - -Kg'^'>) - 


^{X-n)+a 


ab 1 

> eb - 

' 1 

J J 

\ 


V“K-(X -6“) 

2 ^ ^ 


n . 


(BIO) 


Projection onto the unit vector I = laC^ perpendicular to the contact line yields the local force density: 


lar = \-{Kl-Kl)- 


1 / : 9 sin^ ip 


^{X-n)+a 


l- 




n 


2\* - 


-{X ■n)+d- 




•• thcosth cos■)/;sin■)/; -P.,, 

V’+-^-72- 

P 2 


n 


(Bll) 


where the second line only holds for an axisymmetric surface. At the circular contact line: laf^ = — (^ + d) n. 

For the conformally transformed surface one obtains (see Refs. [481450] for further examples of how stresses change 
under conformal transformations): 


laf 


(P \ /e,(l + A2)-2A(A.e,) + 2Axe^)\ y r-n(l - A^) - 2A(1 + n • A) 

I-- 


(B12) 


To analyze the balance of forces let us consider a conformal transformation along the 2 : axis: A = KzZ. The 
projections of the normal and tangent vectors onto the basis vectors are given by 


p ■ Bs = — 2 : • n = cos pj 
p n = z ' Bs = . 


and thus, with n + A = sin^/^p + (A;^ — cos'0)2: one obtains the projections 


P-iaT 


Z ■ iaf 



f COS t/)(l + A^) - 2kz \ 

\ 1 — 2Az cos 0 + A^ y 


-0 sin0 , and 


/ sinV>(l-A^) \ 

\ 1 — 2Kz cos0 + A^ y (1 — A^) 


[cos0(l + A^) - 2 A; 2 ] . 


(B13) 

(B14) 


(B15) 

(B16) 


The axisymmetric solutions we consider are symmetric with respect to the xy plane as well. The two circular contact 
lines at 0+ = 90° + 7 and 0_ = 90° — 7 are transformed into two circles of different length. To calculate the total 
external force let us determine the new positions and lengths of the contact lines. Before the transformation we have 
X = (sin0,O, — cos0)^. After the transformation: 

X = \ sin 0 p + [A(A ;2 — cos 0) + A^] 2 : , (B17) 

where A := (1 — A^)/(l — 2KzCo^p + A^). One obtains for the new contact angles: 0 = 90° + 

arctan [A(A ;2 — cos0) + A; 2 )/(A sin0)]. The lengths of the contact lines are thus: Z = 27rsin0 = 27rAsin0. 

The horizontal components of the force will equilibrate each other due to the axisymmetry. The vanishing of the 
force in the vertical direction yields a condition which can be used to find Kz as a function of the known parameters 
P, d, and '0: 

0 = 2 : • (27rA+sin0+(/0/'')++ 27rA_sin0_(/a/'')_) 

/ -Ccos7(l - A^)^ \ [- sin7(1 + A^) - 2A^] / -Ccos7(l - A^)^ \ [sin7(l + A^) - 2A^] 

\(1 + 2 A ;2 sinq + A2)2y (1 + 2 A ;2 sinq + A^) \ (1 — 2A^ sinq + A^)^ y (1 — 2 A ;2 sinq + A^) 


<^0 = 
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where we have defined C := 
simplifies to 



For 7 = 0, i.e., the solutions which are only confined in one circle, this condition 


0 = ^(1 - A2)2+ 2A^(1 + A^) . (B18) 

A similar equation can be obtained for A = K^x. We observe that the total force acting on the membrane due to 
the container is only equilibrated for special values of A. This implies that the conformal transformation defined in 
Eq. ^ does not yield an equilibrated surface in general. However, we have shown in Sec. |v| that these surfaces can 
be used as initial conditions for subsequent finite element simulations in order to obtain equilibrium shapes. 


Appendix C: Finite element simulations 

In the simulations the fluid membrane is discretized by a set of triangles with N ^ 4000 nodes using subdivision 
finite elements [SUEg. The position of the surface in is parametrized with the local coordinates si and 52 and 
interpolated by 


N 

Xh{s^,s'^) = , 

a=l 


(Cl) 


where Xa is the position of node a and are the Loop subdivision trial functions [53]. We search for equilibrium 
solutions of the energy functional 


E = 


/ 


2H^ + 


Ha 

2 


(V« - Va) 


2 


ds^ds^ + ^iV-Vf , 


(C2) 


where the constants are set to jiA = 10^ and /iy = 5 • 10^ to ensure that the membrane adopts the prescribed target 
values for area and volume to a numerical error of about 10“^. Note that the area constraint is recast in a local form 
{local area constraint method)^ in which ^/a is the surface Jacobian such that A = f ^/ads^ds^ [CT [52] . 

The energy can be expressed in terms of the nodal positions Xa to determine the gradient with respect to x^. This 
yields the corresponding nodal forces due to the elasticity of the membrane (for more details see appendix B of 
Ref. (El). The contact between the membrane and the spherical container is modeled by a harmonic force which is 
only applied to the nodes which leave the container. The force acting on such a node a is given by 

fa = kdln , (C3) 

where k = 1.5 • 10^ is the stiffness constant, da the penetration depth, and n the inward-pointing normal of the 
container. The total force at node a is the sum of all contributions 


f ^ fM 

J a J a ' J a 


(C4) 


To And an equilibrium solution, we integrate these nodal forces in time according to Newton’s equations of motion 
until a balance of forces is reached, z.e., until /^ = 0 (for more details we refer again to appendix B of Ref. [T5]b 
All our input meshes are obtained from the axisymmetric Runge-Kutta solutions (see appendix |^. Every cross 
section of such a solution, which contains the symmetry axis, yields the same profile regardless of the azimuthal angle 
<p. In order to obtain a three-dimensional mesh we discretize this profile in a way which does not depend on <p. We 
first choose the angular stepsize A(/), z.e., the angle between two adjacent planes. Then we establish a meridional 
stepsize, z.e., the distance between two successive points in the same plane. Eor the first point, we calculate the 
distance between this point and the corresponding first point in the adjacent plane di ^2 = This distance is 

taken as the stepsize between the first and the second point of the same plane. We repeat this operation for all the 
points i, di^i-^i = piAcj) until we reach the first point again to close the profile. We obtain a mesh composed of squares 
with edges that are the bonds between two successive points of the same plane and the bonds between two analogous 
points (of same arc length) in two adjacent planes. Dividing these squares into two we obtain triangles that are almost 
equilateral. Moreover, they have the advantageous property that they are smaller close to the symmetry axis where 
the curvature becomes larger. This property is important in order to increase the precision without increasing the 
number of the nodes dramatically. Eor the non-axisymmetric meshes we apply a conformal transformation to each 
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node of the axisymmetric mesh while keeping the same connectivities. In this way we preserve the proportionality 
between the curvature and the size of the triangles. 
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